Optimum Phase Space Probabilities From Quantum Tomography 
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We determine a positive normalised phase space probability distribution P with minimum mean 
square fractional deviation from the Wigner distribution W .The minimum deviation, an invariant 
under phase space rotations, is a quantitative measure of the quantumness of the state. The positive 
distribution closest to W will be useful in quantum mechanics and in time frequency analysis . 
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1. Quasi-probability distributions in Quantum 
Mechanics and Time Frequency Analysis. 

The Wigner quasi-probability distribution W [l]], first 
proposed to calculate quantum corrections to thermo- 
dynamic equilibrium, is now widely used in quantum 
mechanics, statistical mechanics, and technological areas 
such as time-frequency analysis of signals in electrical 
engineering and seismology [2J. The W distribution and 
other quasi-probability distributions such as the Husimi 
Q function [3[ ,the Glauber-Sudarshan P function and 
their s-parametrized generalizations [J can be obtained 
in quantum optics by measuring probability distributions 
of quadrature phases and makingan inverse Radon trans- 
form, i.e. quantum tomography [5j. The Wigner function 
has the unique distinction of being the quantum ana- 
logue of the classical Liouville phase space distribution 
since its marginals reproduce quantum probability den- 
sities of position coordinates %, momentum coordinates 
Pi and indeed of quadrature phases q^ cos 6i + pi sin 8{ for 
all 8i with i taking N values for a 2iV-dimcnsional phase 
space. In time frequency analysis too W has the correct 
marginals reproducing energy densities in time or fre- 
quency. 

Unlike the classical Liouville density, W cannot be in- 
terpreted as a joint probability density, because there 
are quantum states for which W is not positive definite. 
Similarly in time-frequency analysis W, has marginals re- 
producing the energy densities in time or frequency but 
cannot be interpreted as their joint density ; for that 
one uses the positive definite 'Spectrogram' even though 
it does not have the correct marginals. In quantum me- 
chanics, the main reason for the importance of the Husimi 
function Q ( a smeared W function) is that it is positive 
definite ; secondlyas shown by Braunstcin, Caves and 
Milburn, it is the optimum of the distributions obtained 
in the Von- Neumann- Arthurs-Kelly model for joint mea- 
surement of position and momentum [3j. However, its 
marginals differ from the corresponding quantum prob- 
ability densities, even when the W function is positive 
definite. This motivates the variational problem seek- 
ing the best possible positive distribution. The positive 
joint probability we find has immediate utility for quan- 
tum mechanics (especially quantum optics) and in time- 



frequency analysis with obvious transcriptions of the vari- 
ables (q,p going to i, w ) as improvement over the Husimi 
Q function and the Spectrogram Pgp(t,uj) respectively. 
2. Positive joint probability distribution clos- 
est to the Wigner distribution and a measure of 
quantumness. Suppose we know W through quantum 
tomography. We seek a criterion invariant under phase 
space rotations to define the positive definite phase space 
probability density 'closest' to the W function and with 
total phase space integral unity, as necessary for a prob- 
ability interpretation. The criterion of 'closeness' must 
be such that it gives back the W function when that is 
positive definite. In 2N dimensional phase space, with 
units H = c = 1, the Wigner function is given in terms of 
the density operator p , 
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where time dependence of the density operator and the 
Wigner function have been suppressed, q op ,p op denote 
the position and momentum operators and the last equa- 
tion facilitates discussion of rotation properties in phase 
space. In quantum optics, 



q op = (o + a))/V2, p op = -i{a - a<)/V2. (2) 
We vary P(q,p) so as to minimise, 

2 _ Jdqfdp(P(q,p) - W{q,p)f 



JdqfdpW(q,p)i 



(3) 



(which is just the mean of the square of the fractional 
deviation (P — W)/W with the weight function W 2 ), 
subject to the constraints, 



JdqjdpP(q,p) = 1; P(q,p) > 0. 



(4) 



We use Lagrange's method of undetermined multipli- 
ers modified to incorporate inequality constraints. The 
above normalization constraint is equivalent to 



dq dp(P(q,p)-W(q,p))=0, 



(5) 

and the expression for a 2 , using Moyal's well known re- 
sult for phase space integral of W 2 [l[ simplifies to 

a 2 = (2n) N J dqj dp(P(q,p) - W(q,p)) 2 . (6) 

This leads to the Lagrangian, 

L = J dqj dp (P(£ p) - W(q, p)) 2 

+2cJdqjdp(P(q,p)-W(q,p)), (7) 

where c is the Lagrange multiplier. Following a method 
used widely by Martin to incorporate inequality con- 
straints |6j, we prove by direct subtraction that a 2 has 
a global minimum when we choose P(q,p) = P m in(,q,p), 
where, 



Pmm(q,p) = Po(q,p) 8(Po{q,p)), 



(8) 



where 0(x) is the Heaviside 9 function, being unity when 
the argument is positive and zero otherwise, and 



Po(q,p) = W(q,p)-c. 



(9) 



Denoting by L and L m i„ respectively the values of the 
Lagrangian for an arbitrary P(q,p) satisfying the con- 
straints, and by P m in{q,P) , we obtain, 



L - L mm = I (P - P ) 2 dqdp 

JPa>0 

(P 2 - 2PP )dqdp> 0, 



I 



(10) 



Po<0 



since each of the two integrands is non-negative. We com- 
plete the proof by showing the existence and uniqueness 
of a constant c satisfying the normalization constraint, 



W(q,p)-c>0 



(W(q,p) — c)dqdp= 1. 



(11) 



First, if W is non- negative, c = is the unique solu- 
tion, and gives a 2 = 0. Suppose now that W is negative 
in some regions of phase space. The left-hand side in- 
tegral is then > 1 for c < 0, decreases monotonically 
as c increases to positive values until it equals when 
c = max^pW(q,p). Hence there is a unique solution for 
c in the interval [0, max^W(q,p)}. Using this value of c 
we compute the optimum phase space probability distri- 
bution as well as the minimum value of a 2 , an index of 
quantumncss of the state . 



3. Incorporating additional rotationally invari- 
ant constraints in phase space . In addition to the 
phase space volume, the surface of the sphere with centre 

cfchPcl 

(q-qd) 2 + {p-pd) 2 = x 

is an invariant under rotations in phase space ; the ex- 
pectation value of the left-hand side given by the Wigncr 
function agrees with the quantum result and hence may 
be used as an additional constraint. Further, if W re- 
mains positive in the region x > x ma x, we may choose 
P(q>P) = W{q,p) in that region, and for sufficiently large 
Xmax , still find a solution P(q,p) that minimises a 2 under 
the positivity constraint P(q,p) > 0, the normalisation 
constraint, 



//,. 



dqdp(P(q,p)~W(q,p))=Q, 



(12) 



and the additional constraint, 



dqdp(P(q,p)~W(q,p))x = Q. (13) 



If we choose q c i,p c i as the quantum expectation values of 
qop,Pop, the last equation imposes the sum of quantum 
dispersions (Aq) 2 + (Ap) 2 as a constraint. We then prove 
as before that the solution minimising a 2 is, for x < x max 



where 



P mt m(q,p) = Poi(q,p) e(P 01 (q,p)), (14) 



Poi (q, p) = W(q, p)-c-xd, (15) 



provided that constants c, d are found satisfying the two 
equality constraints given above. 

4. Optimum positive joint probability distri- 
butions and Husimi distribution for generalized 
coherent states . The Husimi Q function in 2N- 
dimcnsional phase space is, 

Q{q,p) = (2irr N (a\p\a) (16) 

where |a) are the coherent states, 

a\a) = a\a) , a = (q + ip)/v2., , (17) 

Generalized coherent states Q are displaced excited eigen 
state solutions of the time dependent Schrodinger equa- 
tion for the one dimensional oscillator whose probability 
density packets move classically with shape unchanged, 
and have uncertainty product AqAp = n + 1/2, 

(q\Ht)) =(q- q c i{r)\n) expH(n + l/2)r) 

Cxp(iq cl (T)(q-l/2q cl (T))), (18) 

where, \n) is the n-th excited state and q c i has classical 
motion 



t = cot, q c i(j) = ylcos(r + 



(19) 



The quantum expectation values for position and mo- 
mentum operators are, 



aePllblackPHimmigr 



(<7o P ) = <7d( T ), (Pop) = Qci(t) = Pel- 



(20) 



Wigner functions and Husimi functions can be seen to 
depend on q,p only through the combination, 



x = (q- q c 



(P-Pc 



(21) 



For n = the optimum phase space probability density 
is just the Wigner function which is positive definite. For 
n = 1, 2 the W n (q,p) and Q n (q,p) functions are given by, 



Wi = (2/7r)(a;-l/2)exp(-x), 
Qi = (x/(4tt)) exp (-x/2), 



(22) 



Red-Wigner Function 
Black-Optimum Probability 
Blue-Prob. with Aq Ap 
constraint added 
Green -Husimi Function Q 
n =2 state 




FIG. 2: The same plots as in Fig.l for the n=2 coherent state 
. The optimum and Husimi distributions have a 1 — 0.268084, 
and 0.64429 respectively. 



W 2 = (2/7r)((ar - l) 2 - 1/2) exp (-x), 
Q 2 = (x 2 /(Wir))exp(-x/2). 
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(23) 



We have numerically evaluated the optimum phase 
space probability distributions P m i n ,Pmini- With P m i n , 
for n = l,c = .0105338, AqAp = 1.10753, a 2 = .277049; 
for n = 2,c = .0159546, AqAp = 1.72182, a 2 = .268084. 
With P mi „i ,for n = l,x max = 18, c = .0183706, d = 
-.00141308, a 2 = .287684; for n = 2,x max = 15, c = 
.0423479, d = -.00408366, cr 2 = .322254. With the 
Husimi distribution, for n = 1, AqAp = 2, a 2 = .509259; 
for n = 2, AqAp — 3, ct 2 = .64429. We compared the op- 
timum PmimPmini with W, Q distributions in Figs. 1,2. 
We also compared the corresponding position probabil- 
ity densities in Figs. 3,4. The improvement over the 
Husimi function is obvious qualitatively from the figures, 
and quantitatively from the a 2 values. 




Red-From Wigner Function 
Black-Optimum Probability 
Blue-Prob . with AqAp 
constraint added 
Green-Husimi Function Q 
n=l state 



FIG. 3: For the n=l coherent state , the position probabili- 
ties calculated from the optimum joint probabilities are seen 
to be closer to the true probability (given by the Wigner dis- 
tribution) than the Husimi distribution result. 
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Red-From Wigner Function 
Black-Optimum Probability 
Blu e-Prob . with A q A p 
constraint added 
Green-Husimi Function Q 
n=l state 



FIG. 1: For the n=l coherent state , the optimum phase 
space probability distributions with only normalization con- 
straint, and including additional constraints fixing AqAp are 
compared with the Wigner and Husimi distributions as a func- 
tion of x — (q — q c i) 2 + (p — Pel) 2 ■ The optimum and Husimi 
distributions have a 2 = 0.277049, and 0.509259 respectively. 



5. Optimum Positive Phase Space Densities Re- 
producing N + 1 Quantum Marginals . Cohen and 
Zaparovanny [8j constructed the most general positive 



phase space densities reproducing two marginals of W, 
viz. quantum probability densities of q and p . In 2N- 
dimensional phase space, with N > 2, Roy and Singh Q 
noted that in fact N + 1 marginals of W (e.g. for N = 2, 
probability densities of ((71,(72)5 (pi,Q2), (pi,P2)) can be 
reproduced with positive densities; they conjectured that 
no more than N + 1 marginals can be so reproduced for 
arbitrary quantum states, the "N + 1 " marginal theo- 
rem. This was proved later using an extension of Bell 
inequalities [H| to phase space by Auberson et alplf. 
who also derived the most general positive phase space 
density reproducing 7V + 1 marginals; that density is non- 
unique since it contains an abitrarily specifiable phase 
space function. Among the continuous infinity of positive 
phase space densities reproducing N + 1 marginals which 
one is closest to the Wigner Function ? Our method gives 
a straight forward answer; we give the variational answer 
explicitly for N = 2, and indicate briefly the general- 
ization to TV > 2. Find the phase space density P(q,p) 
obeying positivity, minimum mean square fractional de- 
viation from the Wigner distribution , reproducing the 
quantum probability densities of q, and p . Vary P(q,p) 
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Red -From W-Function 
Black-Optimum Probab. 
Blue-Probab, Aq Ap * 
constraint added 
Green-Husimi Function Q 
n=2 state 




FIG. 4: Same plots as in Fig. 3, for the n=2 coherent state of 
the oscillator. 



to minimise the Lagrangian, 

L = f[(P - Wf + (2A(g) + 2/*(p))(P - W)]dqdp, (24) 

subject to the constraints, 

f(P-W)dp = Q 7 f(P-W)dq = 0,P(q,p)>0. (25) 

L is minimised if we choose for P, the function Pq that 
makes L stationary whenever Po is positive, and zero 
otherwise: 



Pmin = PoO{P ), P = W- \{q) - H(p), 



(26) 



where the multipliers A(g),/ti(p) are determined from the 
constraints. As in Sec. 2, we prove by direct subtraction 
that L — L m i n > 0, the only change being the new choice 
of Po = W — \{q) — (i(p)- The constraints yield a pair of 



coupled integral equations to determine X(q),/i(p): 



(A(g) + n(p))dp = - W(q,p)d P , 

Po>0 JP <0 

(\(q) + ll{p))dq = - f W{q,p)dq , (27) 
P o >0 JP <0 



which complete evaluation of the optimum phase space 
dcnsity.For N > 2, the positivity constraint is supple- 
mented by N + 1 marginal constraints, which can, for 
example, be chosen to be the series of probability den- 
sities of (qi, q 2l ..(?„), {pi,q2, ~qn), ~{pi,P2, ~Pn), in which 
each member is obtained by replacing in the previous set 
one co-rdinate by its conjugate momentum. The optimal 
phase space density is again constructed by a Lagrange 
multiplier method which will now involve N + 1 Lagrange 
multiplier functions. 

6. Conclusion. We have proposed a general method 
to find the positive phase space distribution closest to 
the Wigner distribution that can be used in quantum 
optics as well as in time frequency analysis. A mea- 
sure of quantumness emerges. Qualitative and quanti- 
tative improvement with respect to the Husimi function 
is seen explicitly; e.g. for the generalized coherent states, 
the optimum and Husimi distributions have respectively, 
for n = 1 , a 2 = .277049, and 0.509259, for n = 2 , 
a 2 = 0.268084, and 0.64429 . Similar improvements are 
expected in time frequency analysis. In 2A r -dimcnsional 
phase space the optimum positive density reproducing 
N + 1 marginals can be evaluated. 

A preliminary outline of this work was presented at a 
recent lecture [12| • SMR thanks Sumit Das for a remark 
many years ago on invariance under canonical transfor- 
mations, and Aditi Sen De, Ashok Sen and R. Gopaku- 
mar for discussions. 
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